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Abstract 

We derive a direct connection between Monte Carlo time and physical time in terms of physical parameters, using a 
quantum Hamiltonian with a d- dimensional phonon heat bath interacting with a square-lattice Ising ferromagnet. 
Based on the calculated transition rates, we perform dynamic Monte Carlo simulations using absorbing Markov 
chains to measure the lifetimes of the metastable state at low temperatures. We also calculate the lifetimes 
analytically using absorbing Markov chains. The phonon dynamic gives field-dependent prefactors in the lifetimes 
at low temperatures, that are different from the piecewise field-independent prefactors obtained from the Glauber 
dynamic. 
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1. Introduction 

Dynamics of classical Ising spin systems have 
been studied extensively using the Glauber dy- 
namic [1]. The Glauber dynamic was later also de- 
rived by Martin [2] starting from a quantum Hamil- 
tonian consisting of an exchange interaction and a 
coupling to a fermionic thermal heat bath attached 
to each spin. This derivation was performed in a 
weak-coupling limit in which the square of the cou- 
pling constant multiplied by a characteristic time 
becomes an arbitrary finite constant. In this limit, 
the time evolution of the system was proved to be 
Markovian. 



A different dynamic would change the time 
evolution of the system, but not its equilibrium 
properties. Here we study the quantum nearest- 
neighbor Ising ferromagnet with an applied lon- 
gitudinal magnetic field and a linear coupling be- 
tween the Ising spins and a phonon (i.e. bosonic) 
heat bath. We calculated the transition rates from 
one configuration to another resulting from the 
spin-phonon coupling. We then applied these tran- 
sition rates to dynamic Monte Carlo simulations 
and measured the average lifetimes (r) of the 
metastable state [3]. To measure (r), we set the 
initial configuration with all spins up and applied 
a magnetic field H<0. We define (r) as the num- 
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ber of spin-flip attempts until the magnetization 
reaches zero. Exact predictions [4] of the lifetime 
(r) as T^O are given by 

T]h(t)=T(H, J) = %Ji c -2\H\{£ 2 c -i c + l) , (1) 

where the linear critical droplet size is £ c = \2J/\H\~] , 
\x\ denotes the smallest integer not less than x, 
and r is the energy cost of a critical droplet. The 
critical droplet is a cluster of overturned spins 
which is an l c x (£ c — 1) rectangle with one ad- 
ditional overturned spin on one of the long sides 
of the rectangle. This formula is valid for 2J/\H\ 
not an integer and |i7|<4J. Here we show that 
the specific dynamic crucially affects the time 
evolution of the system. 



2. Dynamic Quantum Model 

The total Hamiltonian we use is H = H sp + 
H p h + 7"4p-ph- The spin Hamiltonian H sp and 
phonon Hamiltonian H p h are given by 



H sp = ~J^2 ~ H z X] a i 

(id) 



(2) 
(3) 



where the first summation runs over nearest- 
neighbor sites only. J(>0) is the exchange cou- 
pling constant, crj are the z components of Pauli 
spin operators attached to lattice site j, H z is a 
longitudinal magnetic field, q is the wave vector 
of a phonon mode, ui^ is an angular frequency of 
the phonon mode with q, and ct and Cg- are cre- 
ation and annihilation operators of the phonon 
mode with q. For simplicity, we ignore the di- 
rection of <f, so hereafter we drop the vector 
symbol on q. If we consider only linear coupling 
between spin operators V(a) and strains caused 
by phonons, the general spin-phonon interaction 
for a single spin a at the origin in real space is 

E 9 \j 2NMu lq ( i( l V (°H - iqVH°)c q ), where N is 
the number of unit cells consisting of one spin, 
and M is the mass of the unit cell. When the Ising 
spins are located at positions Rj , the spin-phonon 
interaction Hamiltonian [5] is: 



^sp-ph — i 



(iq o?)( C t-c,)e«^,(4) 



2NMuj q y H JM q 



where A is the coupling strength between the spin 
system and the phonon heat bath. 

With the given spin Hamiltonian, the dynamic 
is determined by the generalized master equation 
[6,7]: 

' ' - -{ P (t) 7 H sp } 

P^nnWm, 



dt h L 



7m'm — 



(5) 



where p(t) is the time dependent density matrix of 
the spin system, to', n, k, and to are eigenstates 
of H sp , p(t) m > m ={m'\p(t)\m), and W km is a tran- 
sition rate from the m-th to the k-th eigenstate. In 
our case, because there are no off-diagonal terms 
in H sp , this generalized master equation becomes 
identical to the master equation which Glauber 
used in his paper [1], but with different W krn . As- 
suming that the correlation time of the heat bath 
is much shorter than the times of interest, we inte- 
grate over all degrees of freedom of the heat bath 
in order to obtain the transition rates. The tran- 
sition rate from the l-th to the k-th eigenstate of 
H S p becomes 



W k ,i = ^-^2\(n q + l,k |W S p_ p h| n q ,l)Y 



q,n q 



x(n q \pp h \n q }5(Ei - E k - ftw q ) 



(6) 



with the energy eigenvalues of H sp Ei>E k . Here 
n q is the average occupation number of the phonon 
mode with q, and p p h is the density matrix of the 
phonon bath. We can calculate the transition rate 
when Ei<Ek similarly. Eventually we obtain 



W k „ 



A 2 



Q vh d+l c d+2 



(Ek-BkY 



(7) 



where d is the dimension of the heat bath, = 2 
(2ir) for d=l,2 (3), rj is the mass density of a unit 
cell, c is the sound velocity, and (3=1 /ksT. The 
two major differences from the Glauber dynamic 
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are the energy term in the numerator and the neg- 
ative sign in the denominator. In the limit T^O, 
the transition rates vanish when Ei=Ek (this can 
occur for \H\—2J, 4 J). The transition rates satisfy 
detailed balance. 



3. Lifetimes of the Metastable State 

We perform dynamic Monte Carlo simulations 
using Absorbing Markov Chains (MCAMC) [3,8] 
at low temperatures with both the Glauber and d- 
dimcnsional phonon dynamics. For 2J< | if |<4 J, 
the critical droplet consists of a single over- 
turned spin, so the n-fold way algorithm [9] (s=l 
MCAMC) gives adequate speed-ups. However, 
when | if | approaches 2 from below or above, the 
s=2 MCAMC algorithm is needed to prevent the 
system from fluctuating between the all-spins-up 
state and the state with all spins up except for 
one overturned spin. For J<\H\<2J, the critical 
droplet has an L-shape formed by three overturned 
spins. This necessitates using s—3 MCAMC to 
prevent the system from fluctuating between the 
state with all spins up except for one overturned 
spin and the state with two overturned nearest 
neighbor spins. Periodic boundary conditions are 
used. For both dynamics, average lifetimes are 
measured over 2000 escapes with the system size 
L=24. The range of temperatures used is between 
T/J=0.04 and T/J={).2. At very low tempera- 
tures for a given field, MPFUN [10] is used for 
high-precision calculations. Since at low tempera- 
tures the lifetime is the inverse of the probability 
that the system escapes from the metastable well, 
the lifetime is written as 



(r)=Ae? r , 



(8) 



where A is a prefactor and T is given in Eq.(l). 
In a given field, the prefactor A and T can be ex- 
trapolated to zero temperature from the measured 
lifetimes : Tln(r) = Tin A + T (Figure 1). 

The exact lifetimes can be analytically obtained 
from Absorbing Markov Chains (AMC) [3,8] in 
principle. For 2J<|ff|<4J, i c =l, two transient 
states arc needed while for J<|ff|<2J, 4=2, 
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Fig. 1. Tln(-r) in units of Monte Carlo steps as a function 
of T at \H\=2J with the (i-dimcnsional phonon dynamic. 
The symbols are our s=2 MCAMC data. The solid lines 
denote extrapolations to zero temperature. The prefactors 
A can be read from the slopes of the lines. 

seven transient states are needed. For fields lower 
than J, a large number of transient states are 
needed so practically it is not possible to compute 
the exact lifetimes from AMC. The Glauber and 
phonon dynamics provide the same energy barrier 
r. However, the prefactor A in the lifetime result- 
ing from the phonon dynamic is different from 
that from the Glauber dynamic. The Glauber 
prefactor is known to be A=5/4 for £ c =l, A=3/8 
for 4=2 [8], and 3/(8(4 - 1)) for 4>2 [11]. The 
d-dimensional phonon dynamic gives a prefactor 
A that has a non-constant derivative with respect 
to | if | and that depends on the dimension of the 
heat bath. For 2J< | if |<4 J, 



A(\H\,d)- 



4(2|if|-4J) d + (8J-2|if|) d 
4(2|if| -4J) d (8J- 2\H\) d 



(9) 



In the limit that d^O, A approaches the Glauber 
prefactor, 5/4. For an infinite-dimensional bath, 
A^oo for 2J<|if|<5J/2 and 7J/2<|if|<4J, 
A->l/4 for |if|=5J/2, A^O for 5J/2<\H \<7J/2, 
and A->1 for \H\=7J/2. For J<|if|<2J, 



A(\H\,d) 



\H\ d + 2{2J -\H\) d 
2d+3\H\ d (2J -\H\) d 



(10) 



In the limit that d — > 0, A approaches the Glauber 
prefactor, 3/8. For an infinite dimensional bath, 
A -► for J<|if|<3J/2, A -► 1/8 for |if|=3J/2, 
and A — > oo for 3J/2<|ff |<2 J. The expressions 
(9) and (10) for the prefactor are exact as T — > 
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Fig. 2. Prefactors A from the d-dimensional phonon and 
Glauber dynamics as functions of \H\ as T — > 0. Filled 
circles (diamonds) denote the prefactors from the phonon 
(Glauber) dynamic when 2J/\H\ is an integer. When 
2J/\H\ is not an integer, the prefactors can be calculated 
using absorbing Markov chains. 



in a laboratory although the actual value of J may 
vary with structure and surrounding environment 
of materials. 



4. Conclusion 

We derived a relationship between Monte Carlo 
time and physical time in terms of material param- 
eters, starting from a quantum Hamiltonian with 
a phonon heat bath. We applied this dynamic to 
the square-lattice Ising ferromagnet and measured 
the lifetime of the metastable state. This dynamic 
gives a very different low-temperature prefactor 
from the Glauber dynamic. 
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0. At nonzero temperatures there are correction 
terms on the order of e~@ s where 5 depends on 
the field. Figure 2 shows the prefactor A vs. \H\ 
for the phonon and Glauber dynamics. The pref- 
actors for integer values of 2J/\H\ are obtained 
from the MCAMC simulation data, while the pref- 
actors in other magnetic fields are from Eqs. (9) 
and (10). We have confirmed that the measured 
prefactors agree with the calculated values for the 
<i-dimensional phonon dynamic. 

The prefactor from the phonon dynamic diverges 
as | H | — >2 J or 4 J because certain spin flips are 
not allowed. For |if|=4J, the probability to flip 
a spin in the all-spin-up state vanishes, while for 
\H |=2 J the probability to flip a spin in the config- 
uration with three nearest-neighbor spins up and 
one down vanishes. At \H\=2J, the energy bar- 
rier is 2(8 J— 2\H |) for the phonon dynamic because 
the system reaches a critical droplet by creating 
two overturned 2nd-neighbor spins or two over- 
turned 3rd-neighbor spins. Figure 1 shows T=8J 
at \H\=2J . However, this behaviour does not oc- 
cur at \H \— J, so the prefactor from the phonon dy- 
namic is finite as \H\ approaches J. For bulk iron, 
the exchange coupling constant J corresponds to 
about 300 T [12]. So \H\=J is too high to achieve 
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